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Abstract 

SARAH is a Mathematics package optimized for the fast, efficient and precise study of supersymmetric models 
beyond the MSSM: a new model can be defined in a short form and all vertices are derived. This allows 
SARAH to create model files for FeynArts/FormCalc, CalcHep/CompHep and WHIZARD/0 'Mega. The newest 
version of SARAH now provides the possibility to create model files in the UFO format which is supported 
by MadGraph 5, MadAnalysis 5. GoSam, and soon by Herwig++. Furthermore, SARAH also calculates the 
mass matrices, RGEs and 1-loop corrections to the mass spectrum. This information is used to write source 
code for SPheno in order to create a precision spectrum generator for the given model. This spectrum- 
generator-generator functionality as well as the output of WHIZARD and CalcHep model files have seen further 
improvement in this version. 



1. Introduction 

Supersymmetry (SUSY) is still one of the best-motivated extensions of the standard model (SM) of 
particle physics. Even the minimal supersymmetric standard model (MSSM) predicts many new particles 
[U, Q and has been studied intensively for decades. However, so far no hint of any of these particles has 
been found at the LHC 0, 0, H, H, 0, HI- O n the other hand, there are indications for a Higgs mass of 
around 125 GeV 0,113 • Even if this mass can be obtained within the MSSM or the constrained MSSM 
(CMSSM), the interest in non-minimal models has grown. Models with an extended Hig gs sector can not 



only enhance the mass of the SM-like Higgs b oson [Tl [ Fl Fi ji Fl [TBI, Fl , [H, [H, [H, , but also lead to 



candidates 32, 33 
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a rich collider phenomenology (2ll 22l 231 24l |31| and possibly to new dark matter 



H HUfllill, [H nifTTherefore, it is very interesting to study the 
phenomcnological aspects of many extensions of the MSSM. However, many computer tools, such as those 
to calculate mass spectra or the dark matter relic density, or to perform collider studies, still support only 
the MSSM out of the box. 

The Mathematica package SARAH was created to close this gap between model building and the usage of 
well-known computer tools to get precise numerical results for new models. The main idea of SARAH is that 
a new supersymmetric model can be easily and quickly implemented by just defining the gauge structure, 
particle content and superpotential, along with the gauge symmetry breaking and the corresponding field 
rotations. SARAH uses this information to derive the mass matrices, tadpole equations, renormalization group 
equations (RGEs) and vertices. Expressions for the 1-loop corrections to the 1- and 2-point functions are 
also calculated by SARAH. Since version 1.0, SARAH [3 li B] . Pi pf has supported the output of the vertices 
to model files which can be used with CalcHep/CompHep 47L481] as well as FeynArts/FormCalc 0, IHoj . 
Support for WHIZARD/0 'Mega [Hi], H3| was added in version 3.0. 

In addition, SARAH has become the first available 'spectrum - generator - generator': all information 
derived by SARAH necessary for a 2-loop RGE evaluation and the calculation of the 1-loop corrected mass 
spectrum can be exported to Fortran code which can then be compiled with SPheno [53, 54 1 to generate a 
fully-fledged spectrum generator for a given model. This interface between SPheno and SARAH was the last 
missing link between model building and phenomenological studies. The SUSY Toolbox [55| combines all 
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the possible outputs of SARAH to create an environment for the study of SUSY models beyond the MSSM, 



consisting of the powerful public tools SPheno, CalcHep, MicrOmegas [56], HiggsBounds [57|,|58|], WHIZARD 
and SSP. 

The newest version SARAH 3.1 continues this development. The output of model files in the UFO 
(Universal FeynRules output) format [59| is now supported. Thus far, it has only been possible to create 



UFO model files in an automatized way by using FeynRules [60|, lUIl |62j, |63|. This format can be read for 
example by MadGraph 5 |64j |. Hence it is now possible to study all models already implemented by SARAH also 
with MadGraph. Therefore, MadGraph has recently also become part of the SUSY Toolbox. In addition, there 
are other improvements in SARAH 3. 1: the outputs for SPheno, CalcHep and WHIZARD have been amended, 
the PDG.IX scheme is supported and the handling of gauge fixing terms has been automatized. Finally, 
control files for the LHPC Spectrum Plotter can be written. In section we present the UFO output of 
SARAH before we show the other improvements of the newest version in section [3] 



2. Interfacing SARAH and MadGraph 

SARAH can write model files in the UFO format to implement new models in MadGraph 5. This format 
is also supported by other tools like GoSam [(35| or MadAnalysis 5 [66j ]. and soon by Herwig++ 67, 68|. In 
general, the model files written by SARAH include the entire flavor and CP structure of the model. For a 
model which is already implemented in SARAH, it is rather easy to create the requested output by typing 
inside Mathematica 

« $SARAH/SARAH . m 
Start ["$M0DEL"] ; 
MakeUFO [Options] ; 

The meaning of these three steps is the following: first, SARAH has to be loaded by using either the relative 
or absolute path $SARAH to the SARAH installation. The second step is to initialize the requested model, for 
instance the MSSM via $M0DEL = MSSM. A list of all models which are delivered with the public version 
of SARAH 3. 1 is given in |Appendix A| or can be printed by SARAH using the command ShowModels. For a 
step-by-step explanation how to implement new models, we refer to Ref. [i^ and [HST ]. 
During the initialization of a model SARAH calculates the entire Lagrangian density for the mass eigenstates 
as well as the mass matrices and tadpole equations. When this is finished, the routine to write the model 
files in the UFO format is called. Optional arguments can be used for this routine 

• Eigenstates -> $NAME: this defines for which set of eigenstates the model files should be written. 
For the models included in SARAH, $NAME can either be GaugeES for the gauge eigenstates or EWSB for 
the mass eigenstates after electroweak symmetry breaking (EWSB). By default EWSB is used. 

• Exclude -> $LIST: to speed up the output and to keep the model files short it is possible to suppress 
the output of generic groups of vertices which are not needed. By default $LIST = {SSSS, GGS, 
GGV} is used to exclude the scalar four point interactions and the ghost terms. The other possibilities 
are SSS (three scalar interactions), SSV (two scalar - one vector boson interactions), SVV (one scalar 

- two vector boson interactions), SSVV (two scalar - two vector boson interactions), FFS (two fermion 

- one scalar interactions), FFV (two fermion - one vector boson interactions), WV (three vector boson 
interactions), and WW (four vector boson interactions). 

The output written by SARAH consists of the files 

• particles .py: contains the particles present in the model 

• parameters .py: contains all parameters present in the model 

• lorentz.py: defines the Lorentz structures needed for the vertices 



• vertices. py: defines the vertices 



2 



• couplings .py: expressions to calculate the couplings 

These files are saved in the directory 

$SARAH/Output/$MODEL/$EIGENSTATES/UFO/ 

This directory contains also additional files which are model independent and can therefore be used with all 
models: function_library.py, object_library .py, __init__.py and write_param_card.py. These files 
were kindly provided by Olivier Mattelaer. 

To use the model files with MadGraph 5, it is sufficient to copy all files to $MADGRAPH/models/$NAME. Here, 
$MADGRAPH is the directory containing the local MadGraph installation and $NAME is a freely-chosen name of 
a new subdirectory. This directory name is used afterwards to load the model in MadGraph via 

import model $NAME 

MadGraph has a list of pre-defined names for the particles of the SM and MSSM which are used by default. 
However, it could be that there are conflicts between these names and the names used by SARAH in an 
extension of the MSSM. For instance, h3 is defined in MadGraph as the pseudoscalar Higgs in the MSSM, 
but SARAH uses it in the NMSSM for the third-heaviest scalar Higgs. In order to avoid such clashes, a model 
can be loaded using only the names defined by the UFO files via 

import model $NAME -modelname 

To use these model files for the calculation of cross sections, it is, of course, necessary to provide the 
numerical values of all masses and parameters. The necessary input can be obtained by using a SPheno 
module created by SARAH for the given model, see also section 13.11 The SLHA spectrum files written by 
this SPheno module can directly be used with MadGraph. In this way, the chain SARAH- SPheno- MadGraph 
provides a direct link from model building to collider phenomenology. 

Validation. To validate the model files written by SARAH in the UFO format, a list of 1 — > 2 decays and 
2^2 processes was calculated using MadGraph 5.1.4.6. The results were compared with those calculated 
with the MSSM model files delivered with MadGraph. For all 1 — > 2 decays there was always an exact 
numerical agreement. The results for the 2 — >• 2 scattering are shown in | Appendix B| It can be seen that 
all observed differences arc smaller than the numerical error. 



3. Other new features and improvements in SARAH 3.1 
3.1. SPheno output 

SARAH can use all derived information about a supersymmetric model to write Fortran source code. This 
code can be compiled with SPheno and modules for new models are created without the need to write any 
line of source code by hand. The features of the SPheno modules written by SARAH are a precise calculation of 
the mass spectrum using 2-loop RGEs, and 1-loop corrections to all masses. In addition, routines are written 
to calculate the decay widths and branching ratios for SUSY and Higgs particles. Also the calculation of a 
set of electroweak precision observables is included. All calculations are done with the most general CP and 
flavor structure. This functionality of SARAH has been checked against existing spectrum calculators for the 
MSSM, see Ref. [55[, and for the NMSSM, a comparison with NMSSM-Tools was done [z3]. However, the 
SPheno output works also for more complicated models like the minimal SUSY B — L model [6§| , left-right 
symmetric models 0, 71 1 or models with inverse seesaw mechanism [72|, [73j]. To output the SPheno code, 



the command 
MakeSPheno [Options] ; 

has to be used after starting the model. The different options are 

• ReadLists -> True/False: if the vertices and RGEs for the given model have already been calculated, 
these results can be used and a new calculation is skipped. The default value is False. 
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• Eigenstates -> $NAME : the name of the eigenstates for which the SPheno code should be written. 
The default value is EWSB. 



InputFile -> $FILENAME: basic properties of the SPheno version such as boundary conditions at the 
GUT scale are defined in a separate file. The default name for that file is SPheno. m. 

TwoLoop -> True/False: defines if the 2-loop RGEs should be calculated and included in the Fortran 
code. The default value is True. 

StandardCompiler -> $C0MPILER: SARAH writes a makefile to compile the generated code together 
with an existing SPheno installation. The standard compiler defined in the makefile is usually gf ortran 
but can be changed by this flag. 

For more information in general about the SPheno output of SARAH, we refer the interested reader to Ref. (Hij . 
The main improvements of the SPheno output are in the calculation of the low energy observables and of 
the decay widths. SPheno modules written by previous versions of SARAH calculate Sp, b — > sj, U — > Ijj, 
li — > 31 j, a e ,ii,T and the electric dipole moments of the charged leptons. This list has been extended now by 

• fx — e conversion in nuclei (Al, Ti, Sr, Sb, Au, Pb) based on the results of [75| 

• t — > IP with a pseudoscalar meson P° (tt , 77, rj') based on the results of [76| 

• Z — > IJj calculated and implemented by Kilian Nickel [72} 

We want to pronounce that similar to the previously implemented calculation of li — > 3lj also for [i — e 
conversion and r into meson decays not only the photonic contributions are included. Those are often 
assumed to be dominant, however, it has been shown recently that especially the Z can be very important 
in extensions of the MSSM and have to be taken into account 78, [Tt], [Z^, but also the Higgs penguins 



can give sizable contributions [79j. Therefore, in all calculations the possible contributions of Z and Higgs 
penguins as well as box diagrams arc taken into account. 

In addition, the calculation of all decays has been improved by performing an RGE evaluation of all 
couplings from the SUSY scale to the mass scale of the decaying particle. Previously, the values of the 
parameters at the SUSY scale were used in all decays. Furthermore, the calculation of the loop-induced 
decays of a Higgs particle into two photons and two gluons includes now also the dominant QCD corrections 
based on the results given in Ref. [80( ■ 

3.2. Gauge fixing terms 

If ghost vertices were to be calculated by previous versions of SARAH, it used to be necessary to define 
the gauge fixing terms in gauge. However, the new version of SARAH derives these terms automatically 
using the calculated kinetic terms of the Lagrangian. To this end, the condition is applied that the mixing 
between scalar particles and vector bosons vanishes. Afterwards, the derived gauge fixing terms are used to 
calculate the ghost interactions. 

Since it can happen in models with an extended gauge sector that several Goldstone bosons are a mixture 
of the same gauge eigenstate, for each massive vector boson, the corresponding Goldstone boson has to be 
defined 

{{ Description -> "Z-Boson", 

Goldstone -> Ah[{l}]», 
{{ Description -> "Z'-Boson", 

Goldstone -> Ah [{2}]}}, 
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3.3. Output for CalcHep and WHIZARD 

CalcHep. CalcHep is able to read the numerical values of the masses and mixing matrices from a SLHA 
spectrum file or to take the values defined in the variable file; as an alternative, the SLHA+ functionality 



8l| added routines to CalcHep to diagonalize mass matrices to obtain the eigenvalues and eigenvectors by 
itself. Therefore, we have implemented in SARAH the option to include the mass matrices of a given value 
in the CalcHep model files and to parametrize all vertices by the rotation matrices calculated internally by 
CalcHep. To use this possibility, the CalcHep output has to be started via 

CHep [SLHAinput->False , CalculateMasses->True] ; 

The first option has to be used to disable the default approach that all masses and rotation matrices are 
taken from a SLHA spectrum file. The second option is used to write the necessary routines to the CalcHep 
files to diagonalize all mass matrices. If this option were also set to False, the masses and rotation matrices 
would be expected to be given in the variables file (varsX.mdl) of CalcHep. 

Another, small improvement in the CalcHep output is that SARAH now adds the electric charge of all particles 
in the Aux column of the particle list. This is especially helpful for models with an extended gauge sector 
because CalcHep might not be able to derive this information on its own. 

WHIZARD. The speed of the WHIZARD output has been improved significantly: to write the model files for the 
MSSM, the running time has been reduced by more than a factor of 2 (from 12.2 to 5.0 minutes j]]. For a 
more involved model like the B — L-SSM, the speed improvement is even larger (from 6.9 to 1.3 hours). 

3-4- Generalized PDG numbering scheme 

Recently, a generalized PDG numbering scheme (PDG. IX) was proposed which might be more suitable, 
especially for extensions of the SM which include many new particles (8^, . The basic idea is that the 
main properties of all particles (SM/Non-SM, spin, CP character, B — L quantum number, electric charge, 
SU(3)c transformation) are encoded in a nine digit number. This provides not only information about the 
particles but fixes also the PDG particle code of new particles and hopefully prevents the use of randomly 
chosen numbers, and the confusion that this can cause. Therefore, it is now possible to define two different 
PDGs for a given particle in the model files in SARAH: 

{{ Description -> "Gluino", 

PDG -> {1000021}, 
PDG. IX ->{211110001> 
}}, 

By default, the entries of PDG are used. To switch to the new scheme, either at the beginning of a SARAH 
session or in the model files, the following statement has to be added: 

UsePDGIX = True; 

3.5. Output for LHPC Spectrum Plotter 

The LHPC spectrum plotter by Ben O'Leary can produce plots of the SUSY mass spectrum based on 
the information given in a SLHA output file. The repository is available at 

https : //github . com/benoleary/LesHouchesParserClasses_CPP 



Time taken on a Lenovo Thinkpad X220 with 2.53 GHz 
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For the output it is necessary to provide a second control file in addition to the SLHA spectrum file. The 
control file includes information about the paths to the necessary shell tools (gnuplot, latex, dvips, ps2eps, 
rm, mv) and the DTjtX name associated with a PDG number. In addition, the color and column used for the 
different particles are defined in that file. SARAH can provide such a file which works nicely together with 
the spectrum file written by a SPheno module also created by SARAH. By default it assumes the standard 
paths under Linux, while the color and column of each particle can be defined in particles .m using the 
new option LHPC. For instance, to put the gluino in the fourth column and to use purple for the lines, the 
entry reads 

{{ Description -> "Gluino", 

LaTeX -> "\\tilde{g>" , 
LHPC -> {4, "purple"}, 
... », 

As name for the colors all available colors in gnuplot can be used. The control file for a given set of 
eigenstates of the initialized model is written via 

MakeLHPCstyle [$EIGENSTATES] ; 

and saved in the directory 

$SARAH/Output/$MODEL/$EIGENSTATES/LHPC/ 

It is used together with a spectrum file to create the figure by the shell command 
. /LhpcSpectrumPlotter . exe SPheno . spc . $M0DEL LHPC_$MDDEL_Control . txt 

3.6. Calculation of RGEs 

The calculation of 2-loop RGEs for models with many new interactions can be very time-consuming. 
However, often one is only interested in the dominant effects of the new contributions at the 1-loop level. 
Therefore, a new option to the command CalcRGEs has been added to ignore specific parameters in the 
calculation of the 2-loop RGEs 

CalcRGEs [IgnoreAt2Loop -> $LIST] ; 

Here, $LIST contains the parameters (superpotential, soft SUSY-breaking or gauge couplings) which should 
not be taken into account at 2-loop. 

4. Conclusion 

We have presented the improvements in the new version of the Mathematica package SARAH. SARAH was 
created to simplify the study of SUSY models beyond the MSSM and supported in earlier version the output 
of model files for FeynArts/FormCalc, CalcHep/CompHep, and WHIZARD/O'Mega as well as the source code 
output for SPheno. In addition, the new version of SARAH includes also model files in the UFO format which 
can be read by MadGraph 5. This significantly extends the list of SUSY models which can be studied with 
MadGraph at the moment. In this context, MadGraph 5 has also been added to the SUSY Toolbox. 
We have also shown that the handling of gauge fixing terms has been significantly simplified in SARAH and 
the generalized PDG numbering scheme PDG. IX can now be used optionally. Furthermore, the speed of 
the WHIZARD output has been improved and the CalcHep output includes the possibility to diagonalize the 
mass matrices internally. The calculation of the decay widths and branching ratios in the SPheno output 
has been improved and the list of calculated low energy constraints has been extended. Finally, SARAH can 
write control files for the LHPC Spectrum Plotter. 
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Appendix A. Models 

In SARAH 3 . 1 the following models are included: 

• Minimal supersymmetric standard model: 

- With general flavor and CP structure (MSSM) 

- Without flavor violation (MSSM/NoFV) 

- With explicit CP violation in the Higgs sector (MSSM/CPV) 

- In SCKM basis (MSSM/CKM) 

• Singlet extensions: 

- Ncxt-to-minimal supersymmetric standard model (NMSSM, NMSSM/NoFV, NMSSM/CPV, NMSSM/CKM) 

- near-to-minimal supersymmetric standard model (near-MSSM) 

- General singlet extended, supersymmetric standard model (SMSSM) 

• Triplet extensions 

- Triplet extended MSSM (TMSSM) 

- Triplet extended NMSSM (TNMSSM) 

• Models with i?-parity violation 

- bilinear RpV (MSSM-RpV/Bi) 

- Lcpton number violation (MSSM-RpV/LnV) 

- Only trilinear lepton number violation (MSSM-RpV/TriLnV) 

- Baryon number violation (MSSM-RpV/BnV) 

- jUi/SSM (munuSSM) 

• Additional £7(1 )'s 

- [/(l)-extended MSSM (UMSSM) 

- secluded MSSM (secluded-MSSM) 

- minimal B - L model (B-L-SSM) 

- minimal singlct-extcndcd B — L model (N-B-L-SSM) 

• SUSY-scale seesaw extensions 

- inverse seesaw (inverse-Seesaw) 

- linear seesaw (LinSeesaw) 

- singlet-extended inverse seesaw (inverse-Seesaw-NMSSM) 

- inverse seesaw with B — L gauge group (B-L-SSM-IS) 

- minimal U(1)r x U(1)b-l model with inverse seesaw (BLRinvSeesaw) 
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• High-scale extensions 



— Seesaw 1-3 (SU(5) version), (Seesawl,Seesaw2,Seesaw3) 

- Left/right model (f2LR) (Omega) 

• Non-SUSY models: 

- SM (SM, SM/CKM) 

— inert doublet model (inert) 

Appendix B. Validation of UFO output 

Cross section and numerical error calculated with MadGraph 5 using the model files for the MSSM 
included in MadGraph 5.1.4.6 (ctm, Sm) and UFO model files written by SARAH (as, Ss)- The difference 
was calculated as A = (as — o , m)/o's- 
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ee 


— >• 


A°h 


8 


357 


X 


10- 9 


3.46 


X 


10- 4 


8.359 x 10" 9 


3.31 x 


10" 4 


-2.87 x 10" 4 


ee 


-> 


hZ 


6 


223 


X 


10- 5 


3.46 


X 


10- 4 


6.225 x 10" 5 


3.31 x 


10" 4 


-2.89 x 10" 4 


ee 


-»• 


A a H 


4 


804 


X 


10" 5 


3.46 


X 


10- 4 


4.806 x 10" 5 


3.31 x 


10" 4 


-2.91 x 10" 4 


ee 


— ► 


H-H+ 


1 


502 


X 


10- 4 


6.79 


X 


10- 4 


1.503 x 10" 4 


7.65 x 


10" 4 


-3.79 x 10" 4 


ee 


->• 


did\ 


1 


834 


X 


10- 4 


4.18 


X 


10- 4 


1.835 x 10" 4 


5.02 x 


10" 4 


-3.38 x 10" 4 


ee 


->• 


d 2 dl 


1 


834 


X 


10- 4 


4.18 


X 


10- 4 


1.835 x 10" 4 


5.02 x 


10" 4 


-3.38 x 10" 4 


ee 


-> 


W~W+ 


3 


157 


X 


10- 2 


2.16 


X 


10- 3 


3.16 x 10"2 


2.38 x 


10" 3 


-8.93 x 10" 4 
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3d} 



ee — > d3d 3 



■ d 5 d* 5 



ee 
ee 

ee — > d§d\ 
ee — > d^d^ 
ee — > dftd^ 
dd —> e-2,e\ 
dd — > eze\ 
dd — >■ e^e\ 
dd —> e^e\ 
dd — > egee 

dd — >• (S3 §6 

dd — > u\u"[ 
dd — > M3M3 
dd — > U4UI 
dd — > 

dd — > v,6Uq 

dd_->A°h 

dd — > eie* 

dd^xlxl 

dd^xlxl 

dd^xlxl 

dd^xlxl 

dd^xlxl 

dd^xlxl 

dd_->xlxi 

dd_-+xtXi 

dd_^ X2X2 

dd —¥ xtx.2 



dd_ 
dd_ 
dd_ 
dd 



Xi X2 
hZ 
A°H 
H-H- 



dd — >• d\d\ 
dd_^gg 
dd — > U2U2 
dd ->■ d 2 d'2 
dd —> d^d^ 
dd — >• did% 
dd — > d^dl 
dd 
dd — > gg 
dd — > gg 
uu -> gg 

99^ 
d 2 d^ 
dd\ - 

99 -i 



d 6 dl 



99 

-> 99 
d 2 d* 2 



gg ->■ u z u\ 



1.479 x 10" 4 
4.468 x 10" 5 
4.468 x 10" 5 
4.183 x 10" 5 
1.938 x 10" 5 
1.938 x 10" 5 
2.056 x 10" 5 
4.015 x 10" 6 
5.007 x 10" 6 
5.007 x 10" 6 
1.709 x 10" 5 
2.234 x 10" 6 
6.378 x 10" 3 
6.409 x 10" 3 
6.381 x 10" 3 

6.381 x 10" 3 
6.372 x 10" 3 
2.702 x 10" 6 
4.089 x 10" 9 
2.056 x 10" 5 
1.152 x 10" 5 
3.651 x 10" 6 
3.737 x 10" 6 

9.88 x 10" 5 
6.981 x 10" 8 
9.747 x 10" 5 
3.428 x 10" 5 

3.22 x 10" 4 

9.015 x 10" 5 

2.413 x 10" 5 

2.414 x 10 
3.045 x 10 
2.351 x 10 
2.048 x 10 
9.535 x 10" 2 
6.231 x 10"2 
6.378 x 10" 3 
6.376 x 10" 3 
6.388 x 10" 3 
9.539 x 10" 2 

6.382 x 10" 3 
6.382 x 10" 3 

6.231 x 10" 2 
8.496 x 10" 1 

6.232 x 10" 2 
9.729 x 10" 1 

2.872 

9.16 x 10" 1 
8.595 x 10" 3 
8.789 x 10" 3 



4.32 x 10" 4 
6.85 x 10" 4 
6.85 x 10" 4 
7.13 x 10" 4 
3.03 x 10" 4 
3.03 x 10" 4 
4.43 x 10" 4 
5.62 x 10" 4 
4.85 x 10" 4 

4.85 x 10" 4 
4.96 x 10" 4 

3.03 x 10" 4 
8.07 x 10" 4 
8.07 x 10" 4 
8.07 x 10" 4 
8.07 x 10" 4 

8.06 x 10" 4 
6.83 x 10" 4 
3.46 x 10" 4 
4.43 x 10" 4 
7.75 x 10" 4 
3.48 x 10" 4 

6.65 x 10" 4 

7.66 x 10" 4 
3.98 x 10" 4 
1.65 x 10" 4 

7.86 x 10" 4 
6.16 x 10" 4 
8.52 x 10" 4 
5.73 x 10" 4 

5.7 x 10" 4 
3.46 x 10" 4 

3.46 x 10" 4 

4.4 x 10" 4 
3.81 x 10" 4 
7.91 x 10" 4 
8.07 x 10" 4 
8.07 x 10" 4 

8.06 x 10" 4 
3.9 x 10" 4 

8.07 x 10" 4 
8.07 x 10" 4 
7.91 x 10" 4 
3.11 x 10" 4 
7.91 x 10" 4 

3.47 x 10" 4 
2.71 x 10" 3 
3.29 x 10" 4 
7.57 x 10" 4 
7.46 x 10" 4 



1.479 x 10" 4 
4.47 x 10" 5 
4.47 x 10" 5 
4.183 x 10" 5 
1.937 x 10" 5 
1.937 x 10" 5 
2.057 x 10" 5 
4.016 x 10" 6 
5.008 x 10" 6 

5.008 x 10" 6 
1.709 x 10" 5 
2.234 x 10" 6 
6.378 x 10" 3 
6.409 x 10" 3 
6.381 x 10" 3 

6.381 x 10" 3 
6.372 x 10" 3 
2.701 x 10" 6 

4.09 x 10" 9 
2.057 x 10" 5 
1.152 x 10" 5 
3.647 x 10" 6 
3.738 x 10" 6 
9.879 x 10" 5 
6.982 x 10" 8 
9.744 x 10" 5 
3.43 x 10" 5 
3.223 x 10" 4 
9.02 x 10" 5 

2.413 x 10" 5 

2.414 x 10" 5 
3.046 x 10" 5 
2.351 x 10" 5 
2.049 x 10" 5 
9.535 x 10" 2 
6.232 x 10" 2 
6.378 x 10" 3 
6.376 x 10" 3 
6.388 x 10" 3 
9.539 x 10" 2 

6.382 x 10" 3 
6.382 x 10" 3 

6.232 x 10" 2 
8.496 x 10" 1 

6.233 x 10" 2 
9.729 x 10" 1 

2.874 
9.16 x 10" 1 
8.595 x 10" 3 
8.789 x 10" 3 



5.38 x 10" 4 
7.76 x 10" 4 
7.76 x 10" 4 
8.02 x 10" 4 
3.46 x 10" 4 
3.46 x 10" 4 
5.33 x 10" 4 
6.36 x 10" 4 
5.98 x 10" 4 
5.98 x 10" 4 

5.6 x 10" 4 

3.46 x 10" 4 

8.07 x 10" 4 
8.07 x 10" 4 
8.07 x 10" 4 
8.07 x 10" 4 
8.06 x 10" 4 

7.2 x 10" 4 
3.31 x 10" 4 
5.33 x 10" 4 
7.75 x 10" 4 
3.65 x 10" 4 

7.06 x 10" 4 
7.75 x 10" 4 

4.7 x 10" 4 
1.69 x 10" 4 
7.88 x 10" 4 
6.57 x 10" 4 

1.03 x 10" 3 
6.36 x 10" 4 
6.15 x 10" 4 
3.31 x 10" 4 
3.31 x 10" 4 
5.35 x 10" 4 
3.81 x 10" 4 
8.35 x 10" 4 
8.07 x 10" 4 
8.07 x 10" 4 

8.06 x 10" 4 
3.9 x 10" 4 

8.07 x 10" 4 
8.07 x 10" 4 
8.35 x 10" 4 
3.11 x 10" 4 
8.35 x 10" 4 

3.47 x 10" 4 

2.48 x 10" 3 
3.29 x 10" 4 
7.57 x 10" 4 
7.46 x 10" 4 



-4.8 x 10" 4 
-3.63 x 10" 4 
-3.63 x 10" 4 
-1.33 x 10" 4 
2.06 x 10" 4 
2.06 x 10" 4 
-4.04 x 10" 4 
-2.74 x 10" 4 
-2. x 10" 4 
-2. x 10" 4 
-8.78 x 10" 5 
2.24 x 10" 4 
0. 
0. 
0. 
0. 
0. 

4.23 x 10" 4 
-2.69 x 10" 4 
-4.04 x 10" 4 
-6.25 x 10" 4 

1.02 x 10" 3 
-1.89 x 10" 4 

1.14 x 10" 4 
-1.19 x 10" 4 

2.82 x 10" 4 
-5.77 x 10" 4 
-1.15 x 10" 3 
-5.09 x 10" 4 
-4.14 x 10" 6 
-9.94 x 10" 5 
-2.96 x 10" 4 
-2.98 x 10" 4 
-4.05 x 10" 4 
0. 

-1.28 x 10" 4 
0. 
0. 
0. 
0. 
0. 
0. 

-1.28 x 10" 4 
0. 

-1.28 x 10" 4 
0. 

-5.92 x 10" 4 
0. 
0. 
0. 



99 -» bb 
99 -> tt 
tt — >• e\e\ 
tt -)• t^ej 
tt — > £3 eg 
rf — > e4e| 
rf — > ese^ 
rf -> 
rf -5- e 6 e 6 
rf -5- e 3 e 6 
rf — >• uiu^ 

rf — >• 

rf — >• 

rf — > uau% 

TT -)• U 5 l(g 

rf — > ueMg 
rf -> 

rf — > M6'S3 
TT ->■ X1X1 
tt -> X1X2 
tt -)• x?X° 
rf ->• X1X4 
rf -)• X2X2 
rf -»• X3X3 
rf ->• X4X3 
rf -> xlxl_ 
rf -> xfXi 
rf -> X^X^ 

rf X^X2~ 

rf -> XiXt 
tt — > hh 

TT ->• 

rf -)• hH 

TT -)• ftZ 

rf — > /17 

rf -> 

rf -> 

rf -> A°7 

rf -> 

rf -> 

rf — >■ d\d\ 

tt — > did\ 

tt — \ dzd\ 

tt — > did\ 

rf -> dpdj 

rf -)• d e d^ 

rf -)• d 3 

rf — > ^6^3 

xlXi -> M^-W^- 

X?X? -> ^° 



4.104 x 10" 1 
2.033 x 10" 1 
1.508 x 10" 4 
1.508 x 10" 4 
3.975 x 10" 3 
1.352 x 10" 4 
1.352 x 10" 4 
1.782 x 10" 7 
3.821 x 10 
7.907 x 10 
2.713 x 10 
2.713 x 10 
1.687 x 10 
1.787 x 10 
1.787 x 10 
2.036 x 10" 4 
3.948 x 10" 5 
3.948 x 10" 5 
6.334 x 10" 4 
2.324 x 10" 4 
1.231 x 10" 5 
3.063 x 10" 5 
4.241 x 10" 4 
1.008 x 10" 6 
1.853 x 10" 4 
6.011 x 10" 5 
1.018 x 10" 3 
5.731 x 10" 4 
8.021 x 10" 5 
8.023 x 10" 5 

1.36 x 10" 9 
7.418 x 10" 7 
1.084 x 10" 7 
6.472 x 10" 5 

2.55 x 10" 6 
5.283 x 10" 5 
2.661 x 10" 3 
1.957 x 10" 4 
3.108 x 10" 3 
1.496 x 10" 4 
1.834 x 10" 4 
1.834 x 10" 4 
1.479 x 10" 4 
4.468 x 10" 5 
4.468 x 10" 5 
4.183 x 10" 5 
1.939 x 10" 5 
1.939 x 10" 5 
3.247 x 10" 4 
8.048 x 10" 5 
1.439 x 10" 4 



2.39 x 10" 4 
2.84 x 10" 4 
6.78 x 10" 4 
6.78 x 10" 4 
2.94 x 10" 4 
6.84 x 10" 4 

6.84 x 10" 4 
1.46 x 10" 3 
3.9 x 10" 4 

4.31 x 10" 4 
6.1 x 10" 4 

6.1 x 10" 4 
7.23 x 10" 4 

6.85 x 10" 4 
6.85 x 10" 4 
6.48 x 10" 4 

3.02 x 10" 4 

3.02 x 10" 4 

7.42 x 10" 4 
7.51 x 10" 4 

1.17 x 10" 3 

9.43 x 10" 4 

8.03 x 10" 4 

6.18 x 10" 4 
1.94 x 10" 4 
7.15 x 10" 4 
8.57 x 10" 4 
7.09 x 10" 4 
6.56 x 10" 4 
6.5 x 10" 4 
2.09 x 10" 3 
1.81 x 10" 3 
4.38 x 10" 4 
1.56 x 10" 3 

2.17 x 10" 4 
3.68 x 10" 4 
1.34 x 10" 3 

2.19 x 10" 4 
1.42 x 10" 3 
7.36 x 10" 4 

4.18 x 10" 4 
4.18 x 10" 4 

4.32 x 10" 4 
6.85 x 10" 4 
6.85 x 10" 4 
7.13 x 10" 4 
3.03 x 10" 4 
3.03 x 10" 4 
2.77 x 10" 3 
5.17 x 10" 4 
2.85 x 10" 3 



4.104 x 10" 1 

2.033 x 10" 1 
1.508 x 10" 4 
1.508 x 10" 4 
3.977 x 10" 3 
1.352 x 10" 4 
1.352 x 10" 4 

1.786 x 10" 7 
3.819 x 10" 3 
7.907 x 10" 4 
2.714 x 10" 4 
2.714 x 10" 4 
1.687 x 10" 4 

1.787 x 10" 4 
1.787 x 10" 4 
2.036 x 10" 4 
3.948 x 10" 5 
3.948 x 10" 5 

6.34 x 10" 4 
2.322 x 10" 4 
1.232 x 10" 5 
3.059 x 10" 5 
4.232 x 10" 4 
1.007 x 10" 6 
1.853 x 10" 4 

6. x 10" 5 
1.018 x 10" 3 
5.732 x 10" 4 
8.016 x 10" 5 
8.016 x 10" 5 

1.36 x 10" 9 
7.424 x 10" 7 
1.084 x 10" 7 
6.472 x 10" 5 

2.55 x 10" 6 
5.281 x 10" 5 
2.661 x 10" 3 
1.957 x 10" 4 
3.109 x 10" 3 
1.496 x 10" 4 
1.834 x 10" 4 
1.834 x 10" 4 
1.479 x 10" 4 
4.469 x 10" 5 
4.469 x 10" 5 
4.183 x 10" 5 
1.939 x 10" 5 
1.939 x 10" 5 
3.247 x 10" 4 
8.048 x 10" 5 
1.439 x 10" 4 



2.39 x 10" 4 

2.84 x 10" 4 
7.58 x 10" 4 
7.58 x 10" 4 
3.5 x 10" 4 
7.69 x 10" 4 
7.69 x 10" 4 

1.37 x 10" 3 
3.47 x 10" 4 
4.31 x 10" 4 
6.83 x 10" 4 
6.83 x 10" 4 
8.07 x 10" 4 
7.74 x 10" 4 
7.74 x 10" 4 

7.33 x 10" 4 
3.02 x 10" 4 

3.02 x 10" 4 

8.38 x 10" 4 

7.47 x 10" 4 

9.16 x 10" 4 
9.41 x 10" 4 
7.71 x 10" 4 
6.44 x 10" 4 
1.93 x 10" 4 

7.07 x 10" 4 
9.97 x 10" 4 
7.79 x 10" 4 
6.58 x 10" 4 

6.48 x 10" 4 
2.11 x 10" 3 

1.73 x 10" 3 

4.8 x 10" 4 
1.56 x 10" 3 

2.17 x 10" 4 

4.09 x 10" 4 

1.34 x 10" 3 
2.19 x 10" 4 
1.44 x 10" 3 

8.3 x 10" 4 
4.79 x 10" 4 
4.79 x 10" 4 
5.27 x 10" 4 

7.74 x 10" 4 
7.74 x 10" 4 

8.02 x 10" 4 

3.03 x 10" 4 
3.03 x 10" 4 
2.77 x 10" 3 
5.17 x 10" 4 

2.85 x 10" 3 



0. 
0. 

-5.3 x 10- 5 
-5.3 x 10- 5 
-3.97 x 10" 4 
-1.85 x 10" 4 
-1.85 x 10" 4 
-2.17 x 10" 3 
5.53 x 10" 4 
1.85 x 10" 10 
-1.88 x 10" 4 
-1.88 x 10" 4 
-1.78 x 10" 4 
-1.68 x 10" 4 
-1.68 x 10" 4 
-1.72 x 10" 4 
2.71 x 10" 8 
2.71 x 10" 8 
-8.71 x 10" 4 
8.05 x 10" 4 
-8.45 x 10" 4 

1.11 x 10" 3 

2.12 x 10" 3 

1.05 x 10" 3 
-4.46 x 10" 5 

1.73 x 10" 3 

4.6 x 10" 4 
-2.12 x 10" 4 

6.6 x 10" 4 
8.44 x 10" 4 
-2.94 x 10" 4 
-8.63 x 10" 4 
3.23 x 10" 4 
1.55 x 10" 7 
0. 

2.4 x 10" 4 
0. 
0. 

-2.25 x 10" 4 
-6.68 x 10" 5 
-6.54 x 10" 5 
-6.54 x 10" 5 
-1.42 x 10" 4 
-1.57 x 10" 4 
-1.57 x 10" 4 
-1.13 x 10" 4 
2.06 x 10" 8 
2.06 x 10" 8 
3.08 x 10" 10 
0. 

1.39 x 10" 9 
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x?x?- 


+ x?x? 


9.93 x 10" 6 


8.79 x 10" 4 


9.93 x 10" 6 


8.79 x 10" 4 


0. 


x?x?- 




7.99 x 10" 5 


1.25 x 10" 3 


7.99 x 10" 5 


1.25 x lO" 3 


0. 


x?x?- 




2.664 x lO" 2 


1.4 x 10" 3 


2.664 x 10" 2 


1.4 x lO" 3 


0. 


x?x?- 


*XaxI 


8.693 x 10" 4 


1.42 x lO- 3 


8.693 x 10" 4 


1.42 x lO" 3 


0. 


x?x?- 




4.403 x 10" 4 


3.14 x 10" 4 


4.403 x 10" 4 


3.14 x 10" 4 


0. 


x?x?- 


* xixi 


5.931 x 10" 2 


1.55 x lO" 3 


5.931 x 10" 2 


1.55 x lO" 3 


-1.69 x 10" 5 


x?x?- 


■> X.tX2 


1.562 x 10" 2 


1.73 x lO" 3 


1.562 x 10" 2 


1.73 x lO" 3 


0. 


x?x?- 


*xlxl 


5.636 x 10" 5 


9.07 x 10" 4 


5.636 x 10" 5 


9.07 x 10" 4 


0. 
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